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Abstract—This paper takes a fresh look at two classical 
electromagnetic (EM) shielding problems involving 1) a 
time harmonic analysis of a closed, imperfectly 
conducting spherical shell and, 2) a quasi-static analysis 
of a perfectly conducting hollow sphere with an aperture. 
Previous studies of the EM shielding provided by these 
geometries have concentrated on evaluating the E- and 
H-fields at the center of the shield. While the internal H- 
field in the shielded volume of the conducting shell is very 
close to being constant, the same is not true for the E- 
field, where there can be a significant variation in the E- 
field intensity from point to point within the interior. 


For both of these canonical shielding problems, the 
analysis methods are reviewed and are used to determine 
cumulative probability distributions for the fields within 
the shielded volume. 


In implementing the analysis for the conducting shell, 
the use of scaled Hankel functions is described. This is 
used to avoid numerical difficulties in evaluating the 
spherical harmonic solution for lossy medium. 
Additionally, closed-form expressions for the wave 
expansion coefficients in the spherical coordinate system 
are derived. The analysis of the hollow sphere with an 
aperture likewise permits the determination of the E- 
fields anywhere in and around the sphere. 


Index Terms—EM Shielding, Spherical Shield, Sphere 
with Aperture, Spherical Harmonic Expansions, Scaled 
Spherical Harmonics. 


I. INTRODUCTION 


There is a continuing need to understand and describe the 
electromagnetic (EM) field environment inside a protective 
enclosure illuminated by an exterior source. There are 
several EM standards [1]-[3] that provide measurement 
procedures for obtaining shielding effectiveness parameters 
for enclosures and these standards are often used as 
requirements for the design and procurement of equipment 
protected against external EM fields. 

Most of the standards recognize the fact that the EM field 
within a real enclosure will vary with position and 
polarization. Thus, the test procedures usually involve 
making several measurements of the internal field and 


determining a worst-case shielding estimate. However, due 
to time and budget constraints, it is unlikely that sufficient 
measurements are made to develop a robust statistical 
representation of the internal fields. 


It is possible to use a computational model of an enclosure 
to determine the internal EM field and its spatial variability. 
Such models can involve simple canonical shapes like a 
conducting slab, a cavity bounded by two slabs, a cylinder or 
a sphere [4]-[5]. More complicated models of realistic 
enclosures having apertures and conducting penetrations are 
also possible using a finite-difference time domain (FDTD) 
procedure for solving Maxwell’s equations in and around the 
enclosure [6]. 


The simplest, yet somewhat realistic, model for shielding 
is a spherical enclosure. Unlike the infinite cylinder or one or 
more conducting slabs, the sphere has a finite volume, which 
is more typical of a realistic enclosure. Moreover, the EM 
field in the vicinity of the sphere can be described by 
relatively simple mathematical functions that permit a 
numerical computation of the shielding of external fields. 


Using the spherical wave functions defined by Stratton [7], 
Harrison and Papas [8] have developed expressions for the 
E- and H-field at the center of a thin spherical shield due to 
an incident plane wave excitation. Lindell [9] has examined 
Harrison’s solution near the natural resonances of the sphere, 
and Shastry [10] has analyzed a hemisphere being excited by 
a point dipole source. Baum [11] has also examined this 
problem as a special case illustrating the use of the boundary 
connection supermatrix (BCS) of a sphere. While all of these 
investigators have used a modal expansion technique for the 
solution of the sphere shielding problem, Franschetti [12] has 
employed an integral equation approach. In each of these 
references, only the E-fields at the center of the sphere have 
been considered’. 

Real enclosures have openings, so perhaps a uniform 
spherical model is not the best one to use for understanding 
the behavior of the internal fields. References [13]-[15] have 


! This may be due to two reasons. First, at the center only the n = 1 
spherical harmonic is needed, so that evaluating an infinite series summation 
is not required. Second, it is well-known that the quasi-static magnetic field 
is constant within the sphere volume, and using the center as the B-field 
observation point is a good choice. Perhaps it was thought that this location 
would also be suitable for the E-field. However, as will be shown later, the 
E-field in the sphere varies significantly with position, and the E-field at the 
center is significantly lower than the average value of E within the sphere. 
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described analysis procedures for treating a spherical shell 
with a circular aperture. However, their emphasis is on the 
scattered EM field and not the internal field distribution. 
Sancer [16] has described a frequency dependent dual series 
solution for the internal field in a sphere with a hole, and has 
used this model to determine the quasi-static E and H-fields 
only at the center of the spherical cavity. Casey [17] has 
solved the same problem using a quasi-static dual series 
model for the E-field at the center of the sphere. A more 
general frequency dependent solution for the internal fields 
in a sphere has been described by Senior in [18]. This 
solution involves a modal expansion for the E-fields and is 
based on the earlier work of [13]. The applicability of this 
work is in the resonance region of the sphere. 


The work reported in this paper is a re-visitation of two 
classical canonical shielding problems: 1) a thin spherical 
shell made of imperfectly conducting, and possibly magnetic 
material, and 2) a perfectly conducting hollow sphere with an 
aperture. Section II describes the analysis of the penetrable 
sphere by using expansion of spherical vector wave functions 
in the three regions of the sphere: inside and outside the 
sphere and in the wall material. In performing this analysis, 
closed form expressions for the expansion coefficients are 
determined and a method for eliminating numerical overflow 
errors in evaluating various Hankel functions in the wave 
expansions is described. With a computer program 
developed to evaluate the E- and H-fields anywhere within 
the shielded volume of the shell, a Monte Carlo simulation 
has been performed to generate data suitable for describing 
the cumulative probability distributions (CPD) for the 
internal E and H fields. 


In Section III, a quasi-static model useful for computing 
the internal E-field in a sphere with a hole is reviewed. 
Because the dual H-field problem can be solved from the E- 
field solution in this case (see [16]), only the E-field 
shielding is discussed here. A Monte Carlo simulation is also 
performed for this shield, and the corresponding CPDs for 
the E-field are presented. 


Section IV concludes this paper with a review and 
comments. References are provided in Section V. 


II. EM SHIELDING BY A SPHERICAL SHELL 


In this section, the classical solution for shielding of a 
spherical shell of imperfectly conducting material is 
reviewed and the behavior of the internal and external E and 
H-fields is examined. This solution involves a spherical 
harmonic expansion of the fields with expansion coefficients 
determined from suitable boundary conditions on the sphere. 
This analysis is essentially the same as described by Harrison 
and Papas [8] and others. However, [8] only provides the 
expressions for the field expansion coefficients for the 
internal fields and it concentrates on the fields at the center 
of the sphere. In the development here, we will provide 
closed-form expressions for all field coefficients for all 
regions. In addition, in the present development we pay 
special attention to the machine computation of the spherical 
Hankel functions, which are needed in the solution for the 
fields. In particular, we describe a simple modification of the 


classical spherical wave expansion functions that permit an 
accurate evaluation of the E and H-fields inside and outside 
the shell. 


A. Problem Geometry 


The geometry of the spherical shell illuminated by an 
incident plane wave propagating in the z-direction is shown 
in Figure 1. 
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Figure 1. Illustration of a spherical shell illuminated by an 


incident plane wave. 


The incident E-field is in the x-direction, with the H-field 
in the y-direction. The spherical shell is designated as region 
#1, and has outer and inner radii of a and b, respectively. The 
thickness of the shell is denoted by A = a — b. The shell is 
assumed to have the constitutive parameters uj and e; and 


electrical conductivity oz. The material inside and outside 


the shell is assumed to be free space with no conductivity. 
This material is designated as region #2, with parameters 15, 


€ and o» = 0. The spherical coordinate system is described 


by the usual (r, 0, 4) coordinates, as noted in the figure. 


B. Representation of Fields in Spherical Coordinates 


In region #1, the wave propagation constant is 


km X (1-/) for 01.770 €| (1a) 
and in region 72 
k, = Oy) use, - (1b) 


For use later in this paper, the following ratios are defined: 
um and Q=% (Ic) 

k, My 
As developed by Stratton [7] (page 414), the EM field in a 
spherical coordinate system can be expressed as weighted 
a) 


sums of spherical vector wave functions 7; and nO ; 
o 


m,n orn 
where indices n = 0, 1, --- oo, and m = 0, 1, --- n. The symbols 
e and o denote solutions that are even or odd with respect to 
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the x axis, and the index (i) denotes the type of radial 
function used in the expansion. These vector wave functions 
are 








nip, =F" mehr) A SR gd - Q) 
(Ay z cr) Ses mpg 
n =n(n+1)(-1)" a) P"(cosO) 69° mor 
ay 1 ô fy kr 2 (kr i] ESO oos ng (3) 
kr Ok) sin 
M È T kez qr) ] P" (cose) Si. mph 





=) krsin (kr) Cos 


where the radial function 2 (kr) is chosen by the 


parameter i as 


j,(kr) i21 
(rp | dur) i=2 |, (4) 
^ Vr) MP) i-3 


A” (kr) i24 
and k denotes the propagation constant for the specific 
medium in which the spherical waves are propagating. 
The functions P” (x) in Eqs.(2) and (3) are the associated 
Legendre polynomials’, which according to [19, $8.66], are 


defined from the Legendre polynomials P (x) by 


P(x) 2 C)" 0-3 y" E (5) 


From [8] the x-directed incident plane wave shown in 
Figure 1 may be written as a summation of the spherical 
wave functions as 


=E, e ~ Shox s 











2 oy n+l poo yew 6 
o J) aye mo) + AS | (6a) 
E SX jy 2n+1 [ mt”, 35] . (6b) 


ve n(n+1) 
The scattered n reflected") E and H fields for r > a is 


< 2n+1 
E'- EY jy a; me, + jor | (r= a) (7a) 
Da pep im Pia t JEn, | Ca) 





? Tt is unfortunate that there is an inconsistency in the definition of the 
associated Legendre polynomials in the literature. Some references, such as 


Abramowitz [19] include the (-1)” parameter in the definition, as shown in 
Eq.(5). However other authors, including Stratton, omit this factor, with the 
result that the definition of the vector wave functions of Eqs. (3) and (4) can 
vary from text to text. In this paper, we use Abramowitz's definition and 


have modified the wave functions of Stratton by adding the term (-1)” to 
Eqs. (3) and (4). Butler [20] has surveyed a number of widely used texts for 
their usage of this term with the following results. Those authors that use the 


(-1)" term include Abramowitz & Stegun, R. Harrington, D. S. Jones, 


Magnus & Oberhettinger and S. Schelkunoff. Authors that exclude this term 
include J. Van Bladel, J. Stratton, A. Sommerfeld and W. Smythe. 


oo 


2n+1 
E j)” b! im 
AE J) mem n Mein 


2 n=l 





jan, | eza C9 





H= k, 
C) 


Note that the arguments of the radial functions in Eqs.(6) and 
(7) are (ka) and the leading factor for the H-field ky/@up) is 
simply the characteristic wave impedance in medium 2. 


Inside the material of the spherical shell, the EM fields are 
represented in a similar manner, as 














Es panda] "E 
DI S| p, mop, + Ja, P5, 
n=l n(n-1) 
E =E, Sni (b<r<a) 
< n+ 
t» COD —— d, mM HRAD 
2A J) em n Moin Jfa a] 
(8a) 
eC; 2n+l _ . 
2D) |a, mA — iP, m2, | 
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2 J) "mobs e,l,n EN n 255] 


for (b«r <a) (8b) 
with the radial functions arguments being (K ja). 


Inside the sphere void with r < b, the field representation is 











2n+1 = 
-E n c m? + jbtn, (1) (9a) 
x J) n(n + : [a n m, n n "en ] 
jl; oa re em ^ ud 
with arguments of the radial function being (47a). 
The eight coefficients occurring in — Eqs.(7-9) 
a,b’, p, q,, d,, f,, a; and bf are unknowns’ that can 


be determined by the boundary conditions at the interfaces at 
r =a andr = b of the sphere. These boundary conditions are 
that the tangential components of E and H must be 
continuous through the interfaces, and are 
E' +E) =( E° and (E'+E") =(E° 
(E+E) PE), mé (E+E) E) a 
(H+H), -(m), and (H'+H") -(m), 


(10) 
The same form of boundary conditions also holds at r = b. 

In ref. [8] where the EM fields were calculated at the 
center of the sphere, only the n = 1 case was required in the 
summations, At that location, higher order terms of n in the 
field expansions vanish. For the more general case, however, 
other values of n must be considered. By applying the 
boundary conditions independently for each value of n, one 
can obtain a set of eight equations relating the coefficients 
a’, b; --: b; . These equations are 


j,(k,a)+ a’ h® (ka) = p, hP (ka) - d, h? (Ka). (11a) 


n n 


? These parameters are named in the same way as in ref. [8]. 
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flra j, (ka)] +b; [ kah? «| - 


(11b) 
«a [ka A? ka] tf [eni]. 
Ses rj E 
o L, (kra) + EP (a) | ne 
g, A (ka) fh? (ka) 
IE Jj, (ka| +a; | kah? (a) | - 
(11d) 
p, | kai? (a) | 4, [kaht (ka) | 
p, A (kb) d, (kb) =a. j (kb) ^ (1e) 
ela [kbh Qu] + f, [4 0] l= 
(11f) 
b; [b j, kb] 
Q 
Le k b)+ fh (b) |- bj, Ub) aig 
of p, [kbh (b) | +a, [ bnt" «T | - 
(11h) 


a: [kb j, (k,b)] 


where the parameters « and Q have been defined in Eq.(1c). 
Note that the ' symbol in these expressions denotes the 


derivative function, as | xf w] =d | dx(xf(x)). 
C. Scaling of the Hankel Functions and Expansion 
Coefficients 


In using the above equations for determining the 
expansion coefficients, a problem arises in the evaluation of 
the wave functions in region £l. This is due to the 
exponential growth of the Hankel function 5" (Kr) when 


the argument is complex and the frequency becomes large. 
To see this, re£[19, §10.1.16] provides the following 
representations for the spherical Hankel function of the first 


kind AC (kr) 
n (nk) 


nO kr e" j -—n-l kr 
oa um BP re 


In region #1 the propagation constant k; in Eq.(la) is 
complex, and consequently the leading term in Eq.(12) is 


e (ono) )r/2 e" 04401) )r/2 .(13) 


(-2jkr) ae a23 


ur j 





T(kr)s er = e 


4 


This exponential function becomes unbounded as œ, o or r 
— oo, As an example, consider the case of an aluminum shell 
with a radius of 1 meter. As a relatively low frequency of 
1000 Hz, the exponential growth term in Eq.(13) is 


eVa —. 9^. Thus the use of JA" (ka)or h (kb) 


directly in Eq.(11) becomes impossible at high frequencies 
due to numerical round-off errors — even if it's exponential 
growth may be ultimately cancelled out by a small term 
somewhere in a complicated expression. 

In numerical calculations involving (Kk) it is 
convenient to define the scaled Hankel function of the first 
kind as 

(nk)! 


EM i (14) 
kIT(n— k +1) 


Mans j r^ Y (-2j&)* 


k=0 
so that 


h (kr) = T(r) ÁO (kr) (15) 


Similarly, the spherical Hankel function of the second 
kind 4” ( kr) can be written as 


(n+k)! 


(2) -Jr j «n+l ER UA LES , -k 

h (kr) =e (kr) pros pig) di 
2 1 fo 
= Fas hP (kr). 


Equations (15) and (16) use a common scaling function T(kr) 
together with two scaled spherical Hankel functions! 
h (kr) and h® (kr). These scaled functions 


reasonable accuracy over a wide range of values of the 
complex parameter (kr) and they are provided as an option in 
many special function routines. Note that this is not an 
approximation to the Hankel functions, but rather, just a 
factorization of the functions. 


maintain 


To avoid numerical overflow problems in solving the set 
of equations of Egs.(lla-h) we express the boundary 
conditions and Hankel functions in region #1 using the 
scaled spherical Hankel functions and the appropriate scaling 
factors as 


1 
id 


(17a) 


hP (ka) = T(ka) h? (ka); h? (a) = ——— h® (ka) 


hU (k,b) = T (k, b)Á® (kb); KP b) = — sig e» (17b) 


1 

With these scaled Hankel functions for region #1, the 

boundary conditions in Eqs.(11) can be written in a compact 
matrix form as shown in Eq.(18) on the next page. 


. . . ^ . . 
^ [n this paper, any parameter or function with the symbol is designated 
as a scaled quantity. 
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To observe the effects of using the scaled Hankel 
functions and to provide an alternative to a strict numerical 
inversion of Eq.(18) for each value of the summation 
parameter n, one can develop closed-form expressions for 
the expansion coefficients. Reference [8] provides 


expressions for the coefficients a; and bf for the E-fields 


at the center of the sphere and here we generalize their 
results for arbitrary n and for all regions. It should also be 
mentioned that Baum [11] in his Appendix A also provides 
expressions for the expansion coefficients, but not in stand- 
alone terms. He provides ratios of coefficients, relative to the 
coefficients of the incident plane wave. 


Obtaining a symbolic solution is tedious but it can be 
assisted by using a symbolic solver to invert Eq.(18) and 
obtain algebraic expressions for the coefficients. In doing 
this there are two terms that occur in the denominators of the 
various expressions for the coefficients. These are denoted as 
Den, and Den2, and are given in Eqs.(19a and 19b.) 

Note that the functions Den1, and Den2, depend on the 
ratio of scaling functions T(kjay/T(k,b) =e?" —T(KA). 
The bracketed terms A;, A>, B; and B; contain scaled Hankel 
functions f (Ia), hP (a), AP (kb) and AP? (kb), as well 
j, (Sb) and AP (K,a). These 
bracketed terms can be calculated easily at high frequencies, 


as unscaled functions 


with no problems with numerical overflow. The leading 
coefficient T(k,A) still diverges at high frequencies, but 


not as rapidly as T(k,a) or T(k,b). This divergence must 


be offset by a similar decreasing exponential term in the 
expression for the EM fields, as will be discussed in Section 
ILD.3. 


At low frequencies, or for very thin shell walls |k;A| is 
small and the denominator terms Den1, and Den2, can be 
evaluated directly as written in Eqs.(19). However as |&;A| 
becomes large so that IT(KA) »»1/|T(KA) approximate 


expression for these denominator terms can be used. 


The eight expansion coefficients a”, b” 


n? ^n? 


D, ? qn, ? d, ? Sa ? 
a; and b? are given in closed form in Eqs.(20) — (27), which 


are provided in Table 1. For each coefficient the complete 
form and the high frequency approximation are given. Some 
of these expressions have been simplified somewhat through 
the use of appropriate Wronskian relationships between the 
radial functions. Reference [11] also discusses this 
simplification. If Q = 1 (for a non-magnetic shell material), 
additional simplifications are possible using the Wronskian 
relationships. 


Table 1. Closed form expressions for the eight expansion coefficients. 
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- y QA Qb[Ib jb] - 
icc MCI 
Tb) 2 è j, Gb) kb (kb) | 
q T (kb) (= d DEED EIS, uic [xu ec [xbiiauo] -abso iun Q3) 
n ~ Den2, ka 25 fO) : i ' 
K- Ja (k,b)| kib h, (kb) | (e hO (ka)| kah? (ka) | - QR Qa) ka h (k,a)] ) 
_ T(kb) 
u T(kA) a 
È 2 h(kb) [kb j, (sb)] - 
i kal oj (kb) kb fe" (Kb) | 
i (2) : , m . R Y 
1 1 -j h; (kb) [kb j, (k,b)| T (kA) T(kb) (o; (b) kii (b) | -[k,b j, (k,b)] ie 
d Denl, T(kb , ^ 24 
z enl, T(kb) Oj, (D| kbh (b) | " roce. Pe ae ca 
O A? (ka)| Keak (a) | -AP a) ah? ka) | 
o eae 
T(KA)T(kb) " 
(: i ) QA? (kb)[k,b j, (kb)] — 
K| —— 1 
i Ie Jj dob) bh b] 
: iw i 
f K 1 2 QA (k,b)[ kb j, (k,b)] - | (KA)T (kb) [vi ko [koi] -oksan AC) fs n 
n Den2, T (kb) || ka p pO) í , , 
jae) RRR A] [e KP o alia] aita [eekica] | 
rl I 
T(kA)T(b)" 
Q at 
Q 2j 4 "m a JA kb 
d; Den], [m aA kb is Ak, D| LU (k, » - [kb j, (,b)] AO (k, n) | (26) 
IE A (k, o| [k ah” (k, Y - fi (ka) kya h? (a)] J 
ow (zc) 
ae T 1 ka )\ kb 
K 
b; Den2, (z JE "159 [erani] -O[k,b j, kb] Pan) (27) 
B h (k,a)| ka idu] -OP (a) ka naa] ) 

















Notes: T (kA) = 


e^ C 


T(kb)- 


Den, and Den2, are defined in Eqs(19a-b) 
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D. Evaluation of the EM fields in and around the sphere. 


The eight expansion coefficients shown in Table 1 can be 
used with the wave expansions in Eqs.(6) — (9) to determine 
the EM fields in and around the sphere. These calculations 
will be described briefly below. 

1) Outside the sphere (r > a). 

In this region the EM field is the sum of the incident and 
sphere-scattered fields given in Eqs.(6) and (7), respectively. 
For the incident fields only the Bessel function j (k,r) is 


m ! in Eq.(6). 
For the scattered fields, the ied propagating Hankel 
hP (kr) 


and n? 


e,l,n ° 


used in the vector wave functions m, p and nj 


Ln 


function in the vector functions 


(4) 
o,l,n 
The sums for the incident EM field in Eq.(6a-b) can be 
evaluated with no difficulty. Similarly, Eqs.(7a-b) for the 
scattered fields can be summed directly as the expansion 
coefficients a’ and b; are easily evaluated with no overflow 


appears 


m 


problems. 
2) Inside the spherical cavity (r < b) 


In this region the expansion coefficients a. and b; are 


used in Eqs.(9a-b) to evaluate the fields. The vector wave 


functions m^ and no for this calculation all contain 
1 ol! n 


the Bessel aceon J, (kr) , which is well behaved at high 


frequencies and is easy to evaluate. The expansion 
coefficients are given in Eqs.(26) and (27) in Table 1, and 
the high frequency approximations for both are seen to have 
aterm 1/7(k,A) = e "^^ , which is exponentially decreasing 
at high frequencies. This is as expected, since the internal 
field in the sphere should be highly attenuated at high 
frequencies or with thick shield walls. 
3) Inside the spherical shell material (b « r« a) 


For a field observation point inside the spherical shell 
material, the p,, and q, parameters in Eq.(8) are always 


14) 


multiplied by the A?) (kr) terms contained in the m, , and 


—(4 2 2 i 
n : ) vector wave functions. Likewise, parameters d,, and fy 
gmn 


occur together with 4” (kr) functions contained in im? 
o n 


and An 


ol; 


. In these vector functions, we can extract the 


exponential terms of the Hankel functions and use the scaled 
Hankel functions to compute the scaled vector wave 
functions denoted as jn^^ Low n^ 


and 7 Ss from 
Eqs.(2) and (3). For these scaled vector functions the Hankel 
function scaling terms T(k;r) and 1/T(k,r) are used to express 


the vector functions as 


m? -T(kr)m? (28a) 





n? =T(kr)n® (28b) 
(4) l aw 
mo -—— —m, (28c) 
omn T(kr) omn 
iet ew) (28d) 
m,n T(kr) m,n 
Equations (22)-(25) define normalized expansion 


coefficients p, , ĝ,, d DUE together with the normalizing 


factors T(k,b)/T(k,A) and 1/(T(k,b)T(k,A)). Using these 
expressions, along with normalized vectors wave functions 
in Eqs.(28), and introducing an alternate measure of the EM 
field observation point & inside the shell, relative to the inner 


radius, b, as r = b + & we can express the EM fields within 
the shell material as 











n+l ap. ^o (4) 
Ce ue m+ 
ds c JY pon [ 5, o,l;n jå n n | (29a) 
Eerste 9 [2, m9, 7 f, 85. | 
n=l n(n +1) 
n 2n+1 — jk (A+é) in ) 
en m 
TA k, x c J) n+l) [â n — JP. 0] (29b) 
o n 2n+1 _ ^ 
A Me J) ep PRG e Wi, jd, n 272, | 








for (0< € <a) 


The summations in Eqs.(29) are well behaved and easy to 
evaluate. This scaling procedure permits the computation of 
the E-fields over a broad frequency range, which is 
impossible with the conventional wave expansion of Eq.(8). 


E. Numerical Results 


A computer program was developed to evaluate Eqs.(6), 
(7), (9) and (29) using the closed form expressions for the 
expansion coefficients. To check of these calculations it was 
verified that the various boundary conditions on the sphere 
were met and that the overall behavior of the E-fields was 
correct. 


As an additional check, the E-field and H-field transfer 


functions at the center of an aluminum shell (o = 3.54 x 107 
S/m) with radius a = 0.914 m (36 in) and thicknesses A = 
0.794, 1.587 and 3.175 mm (corresponding to 1/32, 1/16 and 
1/8 inches) were computed. These dimensions correspond to 
the same ones used in [8]. For this comparison, transfer 
functions relating the total E and H-field magnitudes to the 
magnitude of the incident E-field are defined as 


|Ecr.6.6) |A(r,0,9)|/ E, 
5 1(Siemen) 





TE(dB) E | and TH (dB) E 
(30) 

In Eq.(30), the argument of the log function for TE is 
unitless. For TH the ratio H/E has the dimension of Siemens, 
so TH is defined relative to 1 Siemen. This latter transfer 
function definition is not ideal, but it has been used to permit 
a direct comparison with the results given in [8]. Figure 2 
presents plots of the E-field and H-field transfer functions, 
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expressed in dB, at the center of the aluminum sphere with 
different shell thicknesses, A. These plots are identical with 
those in ref. [8], and this serves as a partial validation of the 
present computational procedure. Note that usually the H- 
field transfer function is shown relative to the incident H- 
field, not the E-field, and at low frequencies, H/H, — 0 dB. 


The low frequency limit in Figure 2a is —51.53 dB, which is 
exactly 1/Z, expressed in dB relative to 1 Siemen. 
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Figure 2. Plots of the total E-field and H-field transfer functions at 
the center of an aluminum sphere with a = 0.914 m for different 
shell thicknesses, A. (Compare with Figures 2 and 8 of ref. [8].) 


It is interesting to note that the upper frequency response 
of the transfer functions in Figure 2 is 1 MHz. Even for this 
relatively low frequency, the spherical wave function series 
is difficult to evaluate without scaling the Hankel functions. 
To illustrate the robustness of the present scaled solution at 
higher frequencies, Figure 3a presents a plot of the E-field 


transfer functions for the sphere for frequencies up to 1 GHz. 
The responses are reasonable in appearance, and if one 
carefully examines the transfer function for the thinnest 
shell, a hint of small peaks in the curves are seen at 
frequencies above about 200 MHz. 


Figure 3b plots the transfer function for the thinnest shell 
on a linear scale and these peaks are seen more clearly. 
These peaks are due to the internal resonances in the shell. 
Since the shell material is highly conducting, we expect that 
these resonances will occur close to the classical internal 
resonances of a perfectly conducting sphere. According to 
Harrington [21, 86-2], such internal resonances occur at the 
roots of J,(k5b) = 0 for the TE modes, and at the roots of 


[k5b J, (k2b)]' = 0 for the TM modes. 
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Figure 3. Plots of the total E-field transfer function at the center of 
the spherical shell for higher frequencies up to 1 GHz. 
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Table 2 presents the computed roots for the first four TM 
modes, along with the resulting resonant frequencies for the 
sphere with a radius of about 0.914 m. We note that the 
agreement with the frequencies of the resonances in Figure 
3b is very good. Of course, since the conductivity of the shell 
is so high, the internal E-field is extremely small — 1000 to 
2000 dB down from the incident field. 


Table 2. Interior TM resonances for a perfectly conducting shell. 


























No. kb Freq. (GHz) 
1 2.744 0.143 
2 6.117 0.320 
3 9.317 0.487 
4 12.486 0.653 








The primary motive for this study is to understand the 
behavior of the internal E- and H-fields within the shell. 
Using the same spherical shell as in the previous example, 
the transfer functions for the individual components of the E- 
and H-fields have been evaluated along the radial path in the 
direction defined by the angles 0 = 909 and $ = 459. With 
reference to Figure 1, this path is in the x-y plane at z = 0 and 
at an angle of 45° to the x axis. Figure 4 shows these E-fields 
as a function of normalized radial distance r/b from the 
origin to r/b = 2 for a frequency of f= 100 kHz. 

The important thing to observe from this plot is that for 
both the E- and H-fields inside the shell, there are observable 
spatial variations in the E-field components. For the H-field, 
the H, and Hg components are equal in magnitude (on this 
particular radial trajectory), but there is a definite variation of 
the H, g component. For the E-field, the principal component 


is Eg and its spatial variation is significant. Outside the 


sphere both the E- and H-fields approach the incident field at 
distances of r/b = 2. 
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Figure 4. E- and H-field transfer functions for the r, 0 and ¢ field 
components for the spherical shell of inner radius b = 0.914 m and 
shell thickness A = 0.794 mm, at a frequency of 100 kHz. (These E- 
fields are shown as a function of radial distance on a path defined 


by the angles 0 — 90° and ¢= 459.) 


Perhaps a more intuitive measure of the internal shell 
fields is the total E and H-field transfer functions given by 
Eq.(30). Figure 5 illustrates the behavior of the E-field 
transfer functions in the equatorial (z — 0) plane along 
different trajectories defined by the angle ¢. Parameters are 
the same as in the previous example: b = 0.914 m, A = 0.794 
mm, and f= 100 kHz. 


In this figure, it is clear that the total internal H-field in the 
shell is constant for all practical purposes. However, there 
can be a significant variation of the total E-field inside the 
shell — about 60 dB, or a factor of 1000. Furthermore, we see 
that the value of TE for this shield reported in ref. [8] (at the 
center of the shell) is about — 198 dB. This is clearly not a 
representative measure of the overall shielding provided by 
the shell. 


Another way of visualizing the EM field distribution in 
and around the shell it to plot the total field transfer functions 
in contour plots. This is done in Figure 6 for TE and in 
Figure 7 for TH. In these figures we note that for all practical 
purposes the H-field is uniform within the shell. However, 
this is not the case for the E-field transfer function, where the 
TE at the center is considerably smaller than that elsewhere 
in the shell. 


The frequency dependence of the solutions for the E-and 
H-fields of the shell are shown in Figure 8 for a frequency 
range from 100 kHz to 10 MHz. As these are relatively low 
frequencies for the ~ 1 meter sphere size, we do not expect 
much of a change in the shape of the E-field distribution, but 
with a change in the E-field amplitude within the shell. This 
is confirmed by the data in these plots. Clearly, as the 
frequency increases, the internal field strength is reduced, 
due to the inductive shielding provided by the currents 
flowing in the shell [22]. 
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b. Magnetic field 


Figure 5. Plots of the total E- and H-field transfer functions in the 
equatorial (z = 0) plane along different trajectories defined by the 
angle ¢. (Parameters are b = 0.914 m, A = 0.794 mm, and f= 100 

kHz.) 
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Figure 8. Illustration of the frequency dependence of the E- and 
H-field transfer functions, along a radial path defined by the angles 
0 — 90? and ¢= 45°- for the aluminum shell with b = 0.914 m and A 

= 0.794 mm. 


F. Statistical Description of the Internal Fields 


To get a better quantitative description of the behavior of 
the EM fields within the spherical shell, Monte Carlo 
calculations were conducted for a wide range of frequencies 
on the aluminum shell (b = 0.914 m and A = 0.794 mm). In 
these simulations, the total E and H-fields at about 1000 
randomly selected points within the shell were calculated and 
the response amplitudes were binned into histograms 


representing probability density functions. As an example, 
Figure 9 presents the histograms for the E-fields, computed 
at a frequency of 100 MHz. Clearly the H-field is a constant 
within the sphere, but the E-field has a wide range of values. 
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Figure 9. Example of the histogram functions for the internal E 
and H-fields for the aluminum sphere, at a frequency of 100 MHz. 


From the histogram distributions, cumulative probability 
distributions (CPDs) can be calculated. Such distributions for 
the E-field are shown in Figure 10. These distributions 
represent the probability of a randomly selected point in the 
shell having a TE /ess than the value specified on the x-axis. 


As noted in this figure, at very low frequencies (from 100 
Hz to 10 KHz), the E-field CPDs are virtual overlays. At 
higher frequencies, the shielding begins to improve and the 
average values of the shielding becomes larger, with slight 
changes in the shape of the distributions. 
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Figure 10. Cumulative probability distributions for the E-field 
within the aluminum shell (b = 0.914 m and A = 0.794 mm), shown 
for various frequencies. 


As expected, the H-field CPDs are much less interesting, 
due to the almost uniform nature of the internal H-field in the 
shell. Figure 11 presents these results, where there is a slight 
hint of a change in the slope of the CPD at high frequencies. 
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Figure 11. Cumulative probability distributions for the H-field 
within the aluminum shell, shown for various frequencies. 


One way of summarizing the distributions of the internal 
field is by the mean value and standard deviation of the 
histograms. These quantities have been computed for the 
aluminum shell and are listed in Table 3. In addition to the 
data computed for the E-fields within the entire shell volume, 
the shielding values for the E-field and H-field at the center 
of the sphere from [8] are also listed. 

From this summary, it is clear that the use of the H-field at 
the center of the shell as a representative sample of the 
shielding anywhere in the shell is a good measure. However, 
the same is not the case for the E-field at low frequencies, s 
ay below 1 MHz. 
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Table 3.Summary of the mean values and standard deviations of the distributions computed for the 
internal E- and H-fields of the aluminum shell. 


Gace Frequency | Avg. TE | TE Std. TE (dB) 
(Hz) (dB) Dev. (dB) |from Ref.[8] 


1.0E+02 


Avg. TH 
(dB) 


TH Std. 
Dev. (dB) 


TH (dB) from 
Ref.[8] 





1.0E+03 





1.0E+04 





1.0E+05 





1.0E+06 





1.0E+07 





1.0E+08 











III. Low FREQUENCY EM SHIELDING BY A SPHERICAL 
SHELL WITH AN APERTURE 


As mentioned earlier, a completely closed spherical shield 
is not the best geometry to use for trying to represent the 
shielding of a realistic enclosure, because most enclosures 
have apertures or conducting penetrations. To better 
understand the effect of an aperture on the shielding of a 
spherical enclosure, the quasi-static models of Casey [17] 
can be used. These models are useful for electrically small 
enclosures, where k7b < 1. For higher frequencies, the 


dynamic model of ref. [18] can be used, but this is not 
discussed further here. 


Figure 12 shows the geometry of an infinitesimally thin, 
perfectly conducting spherical shell of radius b with an 
aperture, which is immersed in a quasi-static E-field. Two 
different orientations of the E-field are considered: Case 1 is 
with the E-field oriented in the Z direction, and Case 2 is 
with the E-field in the X direction. The aperture is located 
symmetrically at the bottom of the sphere and is defined by 
the opening half-angle 6,, or by the angle a= t — 6,. 


For this sphere with aperture, we are interested in 
computing the internal E- and H-fields and describing their 
statistical distributions throughout the shielded volume. In 
[17], both the E-field excitations and the dual H-field 
excitations of the sphere have been formulated. Reference 
[17] shows that the E- and H-field solutions for this 
geometry are related, with Case 1 for the E-field having the 
same solution as Case 2 for the H-field, and vice versa. In the 
development in the present paper, we will consider only the 
behavior of the E-field for the two cases’. 


$ This is not to underemphasize the need for knowing the internal H-field, 
however. As noted in ref.[23], the excitation of an internal wire in the sphere 
with an aperture requires knowledge of both the E- and H-fields. 











Perfectly 
Conducting 
Shell 





Aperture i d 
‘ ii Case 1 
' f . E 
Quasi-static E-field "S 
p Case2 
par 


Figure 12. Geometry of a thin, perfectly conducting sphere having a 
circular aperture and illuminated by a quasi-static E-field. 


In using this quasi-static shielding model it is important to 
keep in mind that it is inherently different from the 
penetrable spherical shield model discusses in Section 2, 
where at low frequencies the H-field is able to diffuse easily 
into the shield. In the present model with an assumed 
perfectly conducting shield, this diffusion mechanism is not 
present and the shield does provide shielding of the H-field. 


A. Case | -- Axially Symmetric Excitation of the Uncharged 
Sphere with Aperture 


For the case of the z-directed excitation E-field, the 
solution is independent of the $ coordinate. From quasi-static 
considerations, the electric potential function for this 
"incident" E-field is 


14 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


V'(r,0) - —E,rcosO (31) 


and the potential function arising from induced charges on 
the sphere with hole can be expressed as a sum of static 
functions as 


V^(r,0)- Eby, a® B P.(cos@) (r&b) 


n=0 


(32) 
m r —-(n41) 

=E b} a B P.(cos@) (rzb) 
n=0 


where P (cos@) is the Legendre polynomial of order n, and 


d are unknown expansion coefficients. The total potential 


is the sum of the two: 
V(r,0) V'(r,9) -V^ (r,0). (33) 
By using boundary conditions on the sphere that 
V(b,0) — V, a constant (but unknown) for 0 < 0 < a and that 
OV (r,0)/Or is continuous over the opening (r = b, a < 0 < 


71), Casey shows that by using Eqs.(31) and (32) the 
boundary conditions can be put into a dual series equation of 
the form 





Ya P. (cos0)- (0x 0 « a) 

m= (34) 
Y Qna? P (cos@) =0 (a «0 € x) 

n=0 


Using the assumption that the sphere is uncharged, Casey 
further notes that the n = 0 term in the series must vanish and 
he then develops an analytical solution for the coefficients 


i" . (See [17] for more details). The resulting solution is 











a® = l1 ( sin(n-1)&æ x sin(n + 2)a 
ny RE n-l n+2 
VE p) 
j| 2 E M on d E no | sin(n + ve) 
p a+sina n n+l 
The E-field is computed from the total potential as 
E(r,0,6)=-VV(r,0,4) 
Gu tors, Vere 09 
=-— f+ 0 +— 
or r 00 rsin ô 


or in component form, 


E,(r,0)= E,cos0 — 
n-l 
p ) P (cos0) 


E an (for r<b) (37a) 
—(n+2) 
E Laroso[e P(cos@) (for r» b) 


= E,cos0 + 


n-l 
n-l 


Ej(r,0) - —E,sin0 — 


E e y P\(cos@)  (forr«b) 


— (37b) 
--E,sin0— 
(n+2) 
E Xe i P'(cos0) (for r>b) 

n=l 

E,(r,0) =0- (37c) 
In Eq.(37b) the relationship 
dP, (cos) _ P'(cos@) (38) 
d0 E 


is used. 


Casey computes the E-field at the center of the sphere, 
where only the n = 1 term contributes to the sum. In this 
case, he gets an E-field transfer function in closed form as 


E" (0) 





TE? = 


o 


1 1 [sin 0, -sin 20, | 
0, — —sin30, 4 


o 


T z —0, *sinO, 


(39) 





where 0, = zt — at is the opening angle of the aperture in the 
sphere shown in Figure 12. 


B. Case 2 — X-directed E-field Excitation of the Uncharged 
Sphere with Aperture 


The potential function for excitation field in x- direction is 
V'(r,0) » —E,rsinÓ cosó (40) 


and the corresponding potential arising from interaction with 
the sphere with hole is 


V*(r,8,9) = Ebcosé rat 7 ) ese) (r € b) (41) 


n=0 


-(n4l) 
= Ecos aP (7) Pl(cos0) (rzb) 
n=0 


In this expression, the associated Legendre polynomials 
P'(cos@) are used and the expansion coefficients a” are 


different from those found for Case 1. Noting that symmetry 
in this case requires that the total potential on the sphere be 
zero and applying the previously mentioned boundary 
conditions, Casey develops another dual series equation for 
the unknown coefficients for this case: 


Ya? P(cos0)-sinü  (0x0«a) 

n=0 (42) 
Yn +a P. (cos0) 20 (a «0 € zx) 

n=0 


with the solution for q‘” being 


15 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


Tl. : 
PE sin(n -l)a *sinna — 
-1 n-l 


zn(n-1) 


D (43) 


n 





sin(n +1)æ — 7 sin(n+2)a 
n42 


The resulting E-field components inside and outside the 
sphere are calculated from Eq.(36) as 


E,(r,0,9) =E, sin@ cos ¢@— 


- n-l 
E, cosó anz) P'(cos 0) 
n=l 
= E,in0cosó * 


(for r « b) 


—(n+2) 
E, cos $y" aaen Z) P'(cos0) (for r >b) 
n=l 
(44a) 
E,(r,0,ġ) = E, cos 0 cos ġ — 


E, cos > ae E J 


n-l 


= E,cos0 cosó — 


D (cos@) 


forr<b 
20 ( ) 


E, cos pY a? 


n=l 


mD 5 
(5 ) 2g (cos@) (for r» b) 
(44b) 
E,(r,.0,0) =—E, sing + 


p, OS of Q 


" sind ^4 
—-—E,sinó + 


Ya (2) 


n-l 


E y P\(cos@)  (forr«b) 


singy 
" sin 


-(n*2) 
E P'(cos0) (for r» b) 


(44c) 


In Eq.(44b) the derivative of the associated Legendre 
function can be calculated from the relationship 
| (45) 


1 | (n+1)cos6@ P? (cos) 
(n—m+1)P" (cos0) 


n+l 
For this excitation, the E-field transfer function at the 
center of the sphere given by Casey as 


T 2 


s6) = — 


I aa (co - 
sin 








o 


sin 0, ; sin 20, vt sin39, | (46) 


C. Numerical Results 


The 7, Üand ó components E-fields in the vicinity of the 
sphere given by Eqs.(37) and (44) have been calculated 
numerically and a total E-field transfer function TE 
developed as in Eq.(30). As in ref. [18], it was found that a 
large number of terms in the summations were required to 
achieve convergence throughout the sphere. Typically, 150 
terms were used for the present calculations. 

Figure 13 presents a comparison of the calculated TE 
functions at the center of the sphere for Case 1 and Case 2 
with those reported by ref. [17] (in Casey's Figure 2). This 


plot depicts the variation of the shielding as the aperture half- 

angle 6, varies from zero (no aperture) to 180° (no sphere). 
The present analysis overlays exactly on Casey's results, 

and partially validates the numerical implementation of the 


solution. Figure 14 presents the same shielding data as in 
Figure 13, but with TE expressed in dB. 
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Figure 13. Plot of the E-field transfer function at the center of the 
spherical enclosure as a function of the aperture half-angle 8,, for 


Case 1 (E in z-direction) and Case 2 (E in x-direction). 
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Figure 14. Plots of the E-field transfer function (expressed in dB) at 
the center of the sphere, as a function of the aperture angle 6,. 


Since Eqs.(37) and (44) are valid inside and outside the 
sphere, it is possible to develop contour plots of the total E- 
field transfer functions that show the E-field leaking into the 
sphere through the aperture. Figure 15 shows contours of TE 
(in dB) in the vertical z-x plane for the two different 
excitation E-field directions. For these particular plots, the 


aperture half-angle is 6, = 109. 
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Figure 15. Contour plots of TE (in dB) in the vertical z-x plane for 
Case 1 (a) and Case 2 (b) excitations. (Aperture half-angle 6, = 


10°). 


D. Statistical Description of the Internal Fields 


A Monte Carlo simulation with 3000 random internal 
points was conducted and the results for the total E-fields for 
Cases 1 and 2 binned. For the aperture half-angle of 6, = 


109, Figure 16 shows the resulting histograms, which are 
typical of those for the other angles. 
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Figure 16. Example of the histogram functions for the internal E- 
field transfer functions for an aperture half-angle of 0, = 10°. 


Calculations were performed for aperture half-angles of 1, 
2, 5, 10, 20, 30, 45 and 90 degrees and the CPDs of the E- 
field transfer functions computed. Figure 17 presents the 
results for Case 1 and the Case 2 results are in Figure 18. 
Also shown in these figures is a close-up plot of the 
responses near 0 dB, which show that there are some points 
within the r = b sphere that have E-fields larger than the 
excitation field. This occurs at points near the rim of the 
aperture, where the E-field 1s large due to the sharp edge. 


Table 4 summarizes these distributions by the mean values 
of the E-field transfer function and the standard deviations. 
Also shown is the TE value at the center of the sphere. While 
this value is not exactly equal to the mean value of TE, it is 
usually within about 5 dB of the mean. Thus, the value at the 
center does provide a useful measure of the average 
shielding provided by this structure. 
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Figure 17. Cumulative probability distributions for the E-field 
within the sphere with aperture for Case 1, shown for various values 


of the aperture half-angle 6, 
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Figure 18. Cumulative probability distributions for the E-field 


within the sphere with aperture for Case 2, shown for various values 
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Table 3.Summary of the mean values and standard deviations of the distributions computed for the internal E-fields in the sphere with 
aperture. 


Aperture 
Angle (deg) 


TE (dB) 
from [17] 


TE (dB) 
from [17] 





-111.9 


-199.7 





-93.9 


-169.6 





-70.0 


-129.9 





-52.0 


-99.8 





-34.3 


-69.9 





-24.3 


-52.7 





-14.9 


-36.0 











IV. SUMMARY AND COMMENTS 


This paper has examined two canonical shielding 
problems with the goal of gaining a better understanding of 
the EM shielding provided by real shielding enclosures. The 
shield considered here was a spherical shell — one being 
made of finitely conducting material (aluminum) and having 
a finite wall thickness, and the second being a very thin, 
perfectly conducting hollow sphere with an aperture. In the 
study described in this paper, the analysis for the first case 
was frequency-dependent, while the analysis for the second 
was quasi-static. 


The reason for choosing these simple shapes was that the 
calculation of the internal fields could be done 
mathematically through the use of spherical harmonics. This 
provides the possibility of evaluating the E- and H-fields 
anywhere inside or outside the sphere. In this analysis, 
closed form expressions for the expansion coefficients were 
found and these do not appear to be generally available in the 
literature. Moreover, a scaling technique was introduced that 
permits the accurate evaluation of the spherical Hankel 
function terms of the wave functions. This scaling is not an 
approximation to the Hankel functions as obtained by [11] 
and others, but is exact. 


The benefit of this type of solution is that a Monte Carlo 
simulation can be used to develop probability distributions 
for the internal fields, which show the variability of the field 
magnitudes inside the sphere. The difficulty, however, is that 
the solution is in the form of an infinite series of factors, 
which, at times, is difficult to sum. Moreover, there are 
numerical challenges in calculating the required Hankel 
functions of complex argument inside the lossy material due 
to numerical overflow and underflow. 


-2.7 





-10.8 








From the formulation and application of the model 
described in this paper, the following observations regarding 
the shielding of the sphere deserve mention: 


1. Most previous shielding studies of spherical shapes 
have concentrated on the behavior of the E- and H- 
fields at the center of the sphere. For a uniform 
sphere made of finitely conducting material, the use 
of the center as an observation point is appropriate 
for the H-field, as this internal field is seen to be 
nearly constant. However, it is not appropriate for 
the E-field, since there can be large variations of 
this field inside the sphere and the E-field at the 
center is very low. This provides a gross 
overestimate of the amount of E-field shielding 
provided by the sphere. 


2. A better way to describe the internal fields within 
the sphere (and any other enclosure, for that matter) 
is through a cumulative probability distribution that 
represents that variation of the internal fields 


3. For the perfectly conducting spherical shell with an 
aperture, there are also large variations of the E- and 
H-fields within the interior. However, it is noted 
that for this type of shield, the value of the field at 
the center of the sphere is close to the average value 
of the field inside — at least for apertures with 
opening half-angles from 0 to 90 degrees. 
Nevertheless, there can be significant variations of 
the internal fields, with a standard deviation of 5 to 
30 dB being noted the cases considered here. 


4. Finally, while not considered here, but very 
important nevertheless, the presence of conducting 
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penetrations into the interior of the sphere will 
radically change the statistical behavior of the 
internal fields. Since this type of penetration is very 
common, the shielding of most practical enclosures 
will be ultimately determined by this type of 
penetration and not by diffusion and aperture 
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